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Abstract 

Physics goals of the STAR Experiment at RHIC in recent (and future) years 
drive the need to operate the STAR TPC at ever higher luminosities, leading to 
increased ionization levels in the TPC gas. The resulting ionic space charge intro- 
duces field distortions in the detector which impact tracking performance. Further 
complications arise from ionic charge leakage into the main TPC volume from the 
high gain anode region. STAR has implemented corrections for these distortions 
based on measures of luminosity, which we present here. Additionally, we highlight 
a novel approach to applying the corrections on an event-by-event basis applicable 
in conditions of rapidly varying ionization sources. 
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1 Introduction 



The time projection chamber (TPC) used by the STAR experiment at RHIC 
has several potential sources of field distortions [1] . While most of these sources 
are static, the buildup of slow-drifting positively charged ions in the volume 
gas generated from standard operation of the TPC varies with the quantity 
of charged particles traversing the TPC, and thereby both the luminosity of 
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the collider and the multiplicity of charged particles emitted by the collisions. 
The variations in this "space charge" can occur on time scales down to what 
it takes the ions to drift the length of the chamber, which is approximately 
one half second for the STAR TPC. 



2 Space Charge Distortions 



Modeling the distortions due to space charge is a straightforward process be- 
ginning with a postulation of the typical three-dimensional distribution of 
ionization in the TPC. The nearest measure we have of this in STAR is a 
record of the distribution of electron clusters reaching the TPC endcap aver- 
aged over many events using a so-called "zero-bias" trigger (which is random 
with respect to collision times, removing any biases related to the definition of 
a collision). This measure integrates out any drift-direction dependencies, but 
compares well in radial dependence (approximately as inverse radius squared) 
for ^/snn = 200 Au+Au collisions to a simulation using the HIJET event 
generator [2]. The simulation indicates a uniform distribution of charge in the 
drift direction. 

We use the HIJET charge distribution integrated along the distance from the 
endcap to any point in space (representing the effect of continual collision con- 
tributions) in conjunction with the boundary conditions of grounded surfaces 
surrounding the TPC gas volume to solve for the electrical potential due to 
space charge. An analytical solution is not achievable, so we use a numerical 
relaxation to solve for the potential on a grid in two dimensions (with assumed 
azimuthal symmetry) and interpolate. An electric field is obtained from the 
potential and is treated as a perturbation atop the normal drift field. The 
distortions to the measured positions of electron clusters are then calculated 
by integrating the effects of this perturbing field (which depend on operating 
conditions of the chamber) along the path from a point in the TPC to the 
endcap where the clusters are measured [1]. The amplitude of this distortion 
is directly proportional to the quantity of space charge (psc) present. In prac- 
tice, we calibrate the average charge density over the volume of the chamber: 
(PscAo)- 

Because the Lorentz force on the drifting electron clusters is proportional 
to the cross product of the electric and magnetic field (aligned along the 
drift direction in STAR) vectors, the principal distortion of consequence is 
azimuthal, and is plotted in Fig. 1. This distortion has the effect of rotating 
reconstructed tracks in the transverse plane about a point midway along their 
path through the TPC 
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Fig. 1. Simulated shape of the potential due to space charge in the TPC (left) and 
the azimuthal distortions of electron clusters (right) caused by drifting through that 
potential as a function of radius R and drift Z. The cathode is at Z=0, and electron 
clusters drift to the endcaps at high Z. 

3 Space Charge Corrections 

Knowing psc is sufficient to subtract the calculated distortions from mea- 
sured electron cluster positions to obtain their approximate original, undis- 
torted positions. In the absence of direct measures of psc ; a measure of the 
distortion to tracks (fit from distorted clusters) may suffice to indirectly de- 
termine psc- Simulation shows that for any given distorted primary particle 
track, its signed distance of closest approach (sDCA) 1 to the collision ver- 
tex is approximately linearly proportional to space charge, and we can obtain 
QrSk = PscVsDCA™ ck , where C t s ™k depends on the locations of points on the 
track. Each real track can then be used to derive an observed space charge: 

obs _ f-ysim riPA ol>s — n sim . ^DrA obs /^DrA sim 

PSC — '-"track SL,bA trad ~ PSC l S U^A track / SU^A track ) 

To understand the scale of this distortion, it is worthwhile to note that some 
recorded events exhibited beyond 1cm offsets in (sDCA), the mean of their 
sDCA° b a s ck distributions. 2 

A distribution of Psc va l ues from any given collision event will include a 
background from secondaries which naturally do not point to the collision 
vertex, and will be smeared by the intrinsic resolution of the TPC to measure 
sDCA. As seen in the distribution from a single very high multiplicity event 
in Fig. 2, the centroid of a peak formed by primaries provides a means to 
determine p§Q more accurately. 



The sign is determined by the Z-direction of the cross product of the track mo- 
mentum vector at its closest approach and the vector pointing to the collision vertex, 
essentially identifying on which "side" of the vertex the reconstructed track passes. 
2 We use only TPC tracks with at least 25 points, pseudorapidity within ±1, and 
transverse momentum between 0.3-2.0 GeV/c for all sDCA and psc measurements. 
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Fig. 2. Observed space charge density (averaged over the volume of the TPC) 
determined from individual tracks in a single high- multiplicity event. The mean of 
a Gaussian peak (formed from primaries) is fit to extract p^Q for that event. 

To be effective, the value of psc used to correct the distortions must be up- 
dated on time scales shorter than the fluctuations caused by collider operating 
conditions. During the 2000 through 2003 years of operating RHIC, scalers of 
trigger counter rates recorded online (during runs) every 30 seconds served to 
measure these fluctuations sufficiently [3]. Along with a significant luminosity 
increase in 2004, however, these fluctuations were observed in the systematic 
behavior of sDCA distributions on sub-second time scales. 

An event-by-event (E-by-E) method using only pg^ 1 from individually recorded 
events suffers from insufficient statistics to get a good measure in most events. 
To compensate, we can take advantage of the fact that psc fluctuations cannot 
occur on time scales much shorter than the drift time of ions in the TPC. We 
do this by building a running sum of Psc? fr° m each event and previous events 
downweighted appropriately by their age. Because we measure p 1 ^? from events 
which have already been corrected with some value p^Q d , we set the new value 
to be p^™ = pgQ d + Psc - This method is self- correcting in that even if the 
conversion factors C t s ™ k are not perfect, psc will quickly converge to a value 
which brings the sDCA distributions to peak at zero. 

Weaknesses in this technique include events at the start of data files (for which 
there are no previous events), sizable time gaps between some events, and 
series of low multiplicity events for which insufficient statistics are obtained 
within short time scales. The first problem is solved by performing a prepass 
on the first few events in each file to determine a viable initial Psc PasS > which is 
then used in the production pass until the E-by-E method becomes applicable. 
The latter issues are handled by falling back to Psc Pass f° r such events until 
the E-by-E method can again be useful. Backgrounds which introduce charge 
distributions different from the HIJET model can also degrade performance. 

Fig. 3 demonstrates that the fluctuations in psc determined by the E-by-E 
method are not artificial. In two independent but concurrent sets of events, 
similar behaviors can be seen on sub-second time scales, while differences 
illustrate the uncertainty on (pscAo) m the method of about 0.0001 V/cm 2 . 
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Fig. 3. Observed volume-averaged space charge density measured and used in the 
E-by-E method for two selections of independent but concurrent events versus time. 

4 Ion Leakage Around the Gated Grid 

Studying residuals of TPC cluster positions from track fits revealed that an 
additional source of ions is also present in the TPC. A discontinuity in the 
residuals at the gap between the inner and outer readout wire chambers of the 
TPC, evident in Fig. 4, is consistent with incomplete blockage (by the gated 
grid) at this gap of ions created in the high gain region around the anode 
wires. This allows a sheet of ions to flow from this gap across the TPC gas 
volume to the cathode. 

Again, we can model the distortions from this leak around the gated grid in 
the same manner as the space charge, providing a map of cluster position 
corrections whose magnitude is proportional to the amount of leaked charge 
(pieak)- These distortions similarly affect sDCA, and pi ea k was found to scale 
with collision rates in the same manner as pso A calibration was performed 
to find the ratio (D) between pi eak and psc which removed the residual discon- 
tinuities while simultaneously zeroing sDCA in a sample of events. And the 
E-by-E correction was modified to track the two distortions together: 

(Psc + PlSc) = (Pic + &) • (sDCA° b : ck /sDCA^ ck ) , p lcak = D ■ p S c 




Fig. 4. Residuals of TPC tracks over R and Z in a selection of events acquired 
during high luminosity before (left) and after (right) leakage distortion corrections. 
The gap between sector wire chambers is at R ~ 122cm. 
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Fig. 5. Performance of the ionization distortion corrections as measured by the 
distributions of (sDCA) (error bars are the spread (RMS), diamonds the mean) 
versus luminosity (represented by the rate of zero degree calorimeter (ZDC) coinci- 
dences [3]) for y / SNN = 200 GeV AuAu collisions. 

5 Summary 

We have identified and corrected for distortions due to ion charge buildup in 
the STAR TPC. With the onset of significant short time scale fluctuations 
in the sources of the ions which were not monitored with fine time granu- 
larity during data acquisition, we have developed a technique to determine 
and adjust for the fluctuations during reconstruction on an event-by-event 
basis. Performance of the corrections can be assessed by examining the dis- 
tribution of (sDCA) as a function of luminosity, shown in Fig. 5. Here we see 
that the spread in (sDCA) is contained to within approximately 1mm at all 
luminosities, and the mean is kept to within a few hundred microns of zero. 
In 2005, online monitoring with one second granularity was implemented and 
will provide further assessment of the technique's success. 
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